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Abstract 

£SJ , We apply the functional rcnormalization-group (FRG) equation to analyze the 

nature of the QCD critical point beyond the mean-field approximation by taking 
into consideration the fact that the soft mode associated with the QCD critical 
point is a linear combination of fluctuations of the chiral condensate and the 
quark-number density, rather than the pure chiral fluctuations. We first construct 
■ an extended quark-meson model in which a new field ip corresponding to quark- 

number density is introduced to the conventional one composed of the chiral 
fields (er, 7r) and the quarks. The fluctuations of the quark-number density as 
well as the chiral condensate are taken into account by solving the FRG equation 
which contains a and ip as coupled dynamical variables. It is found that the 
| mixing of the two dynamical variables causes a kind of level repulsion between 

the curvature masses, which in turn leads to an expansion of the critical region 
| of the QCD critical point, depending on the coupling constants in the model yet 

to be determined from microscopic theories or hopefully by experiments. 
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§1. Introduction 



>: 

C^- . The phase diagram of Quantum Chrome-dynamics (QCD) has been one of the 

main topics of theoretical high energy nuclear physics. In fact QCD phase diagram is 
qq ■ expected to have a rich structure in a temperature (T) and quark-chemical potential 

(fi) plane P-El JH) At low temperature and low chemical potential, the matter is con- 
sidered to be strongly correlated and a naive perturbation does not work. Therefore 
non-perturbative methods are necessary to analyze the phase structure. 

Lattice Monte-Carlo simulation is a first-principle and powerful approach for 
non-perturbative properties of QCD. At finite quark-chemical potential, however, we 
encounter the notorious sign problem which makes the statistical importance sam- 
pling method difficult. There are various proposals to overcome the sign problem,® 
but the cold and dense region fi ^ T is still inaccessible by lattice simulations. 

On the other hand, approaches based on effective models embodying proper 
symmetries with low-energy effective or collective degrees of freedom are legitimate 
and useful methods to investigate the phase diagram where we need to consider only 
low-lying soft modes with the energy scale of temperature ~ 200 MeV. Many inves- 
tigations of the QCD phase diagram have been done on the basis of effective models, 
including (Polyalov loop extended-)Nambu-Jona-Lasinio ((P)NJL) modeP^^®® 
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or (Polyalov loop extended-) Quark-meson ((P)QM) modelP^'EB'nS Whereas there 
is a strong model dependence of the position of the critical point many chi- 
ral effective model calculations suggest the existence of a QCD critical point at 
finite temperature and quark-chemical potential where the crossover transition in 
the lower-/i region becomes a second order phase transition followed by a first-order 
in the higher-^ region. 

In the vicinity of the QCD critical point, there is a critical region in which suscep- 
tibilities^ or their higher derivatives^ 1 show a critical behavior. For experimental 
searches of the QCD critical point, it is necessary to examine a shape and size of the 
critical region as well as the position of the critical point. Indeed the shape of the 
critical region has been also investigated QSJEHIHJiQID 

In general, there should exist at least one softening mode inherent in the second- 
order nature of the transition around a second order phase boundary. Analyses of the 
dynamical soft mode have been done on the basis of the quantum field theory^ 1 '^ 1 
and the Langevin equation.^) Both analyses suggest that the soft mode around 
the critical point should be the quark-number density fluctuation or phonon mode 
(particle-hole excitation) , in addition to the sigma mode as the amplitude fluctuation 
of the chiral order parameter. This situation is identical to the first-order liquid- 
gas phase transitions of water or the nuclear matter. As we know well in the a-oj 
modelP'''^ 1 the sigma mode couples with the quark-number density mode at finite 
chemical potential. As a result, the linear combination of these two modes would be 
realized as the soft mode around the critical point. 

Recently alternative possibilities of the phase structure are suggested, for ex- 
ample, inhomogeneous chiral condensation phases with the Lifshitz pohrPS'GSJ'EZl 1 
or multiple chiral critical points P^'^'® i n the present work, however, we 

focus on the conventional phase structure with a single first-order phase boundary at 
low temperatures and a single critical point. Since the coupling between the quark- 
number density and the chiral condensate is inevitable at finite chemical potential, 
the analysis in the present work is also appropriate for the other candidates of the 
phase structure. 

Near the second-order phase boundary, fluctuations of the soft mode are en- 
hanced and the mean-field approximation breaks down. We need to include the soft 
mode and its fluctuation effects in order to correctly describe the QCD critical point 
and its critical region. The functional renormalization group (FRG) method^ 'ED 1 
is one of the machineries to evaluate an effective action beyond the mean-field ap- 
proximation, see^'^ for a general introduction. Although there are many analyses 
on the QCD phase diagram using FR(gBMbM)mmmm wne re effective models 
are mostly used, there have been only few works on the QCD critical point which 
include the quark-number density fluctuation as well as chiral modes beyond the 
mean-field approximation. 

In the present work, we start from an effective model including the quark-number 
density fluctuation as well as the sigma and pion modes. We apply the FRG method 
to take into account the fluctuations of the soft mode. A key quantity in our analysis 
is a strength of the coupling between the sigma and quark-number density. We 
thereby examine not only the phase structure around the critical point but also the 
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sensitivity of the critical region to the sigma-density coupling strength. 

This paper is organized as follows. In Sec. [21 an evolution equation for the scale 
dependent effective potential is derived. In Sec. 12.11 we introduce an effective model 
which should describe the dynamics near the QCD critical point. The new field <p is 
introduced in order to describe the quark-number density. The FRG method and a 
suitable approximation for our model is discussed in Sec. 12.21 Results of numerical 
calculations are given in Sec. [3 Section 0] is devoted to the summary and outlook. 

§2. Formulation 

2.1. Effective model near QCD critical point 

In this section we introduce an effective model for describing the two-flavor QCD 
near the critical point. Our starting point is the quark-meson (QM) modeP* which 
is one of the chiral effective models of the low-energy QCD. Since the QM model 
explicitly contains the meson degrees of freedom in the Lagrangian, we can easily 
incorporate the fluctuations of the mesonic modes in the framework of the FRG. The 
Lagrangian of the QM model for 2-flavor and 3-color degrees of freedom reads 

£qm = - 9s(o- + rysr ■ vf)]^ + ^(<V) 2 + \i d ^f ~U(a 2 + tt 2 ) + ca, (1) 

with the mesonic potential U(a 2 + tt 2 ) = ^-(a 2 + tt 2 ) + ^(a 2 + tt 2 ) 2 . The explicit 
symmetry breaking term (ca) corresponds to a finite current quark mass. Except 
the explicit symmetry-breaking term, the Lagrangian is invariant under the chiral 
SU(2) L x SU(2) R transformation. 

The QM model has been extensively applied to describe the global structure of 
the QCD phase diagram, i.e., a chiral symmetry breaking at low temperature and 
low chemical potential and its restoration at high temperature and/or high chemical 
potential. However, near the QCD critical point, the QM model misses a relevant 
degree of freedom: the fluctuation of the quark-number density. In fact, the soft 
mode associate with the QCD critical point is a linear combination of the sigma field 
and the quark-number density J25U2D122J Thus the QM model should be extended to 
correctly describe the thermodynamics near the QCD critical point. 

We first introduce a new bosonic degree of freedom p following the analyses in 
RefsP andP the new field ip corresponds to the quark-number density with an 
appropriate normalization and naturally couples to ^"fotp. Thus our effective model 
takes the following form: 

L = <ij)[i@ -g s (a + ij 5 T ■ tt) + g d <P7o]^ + ^(^<t) 2 + \{d^Tr) 2 + ^-(d^ip) 2 

2 (2) 

m 

Whereas the effective model ([2]) looks similar to the NJL model with a vector 
interaction,H3'^'S3'ESJ'!^'!^ our model has only the mode corresponding to the 
zeroth component of the vector fields. One might suspect that the other compo- 
nents of the vector fields should be included in order to respect the manifest Lorentz 
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invariance. It should be, however, noted that we apply the model to finite chemical 
potential where baryonic matter is formed and the Lorentz symmetry is explicitly 
broken. Moreover the other components of the vector fields are irrelevant to the 
critical phenomena near the QCD critical point because we do not have finite con- 
densations of the spatial components and then the spatial components is decoupled 
from the chiral modes. In this sense, our model embodies enough ingredients for the 
investigation of the critical phenomena associate with the QCD critical point. 

There also exists a difference in the nature between the ip and the zeroth com- 
ponent of the vector fields. In our set up, the field •p which is supposed to describe a 
particle-hole excitation and causes an attractive interaction between two quarks^ 
while the zeroth component of the vector field is usually introduced as a source of 
repulsive interaction. We will come to this point later. 

At finite chemical potential, both the chiral condensate and quark-number den- 
sity take finite expectation values, {a} ^ and (p) ^ 0. We calculate the effective 
potential with the fluctuations of the chiral modes and tp by solving the FRG equa- 
tion and determine the expectation values of these modes from the minimum of the 
effective potential. 

2.2. Functional renormalization group 

In order to evaluate the effective potential with the fluctuations of the quark- 
number density as well as the chiral modes, we adopt the FRG equation. We make 
an approximation suited to solve the QM model and derive an evolution equation 
for a scale dependent effective potential. 

In the FRG approach, we introduce a scale (fc) dependent effective action 
which should satisfy the following properties at ultraviolet (k = A) and infrared 
(k = 0) scale: 

r " =S (3) 

where S and r are the classical action and quantum effective action, respectively. 

We can construct the scale dependent effective action by inserting a mass-like 
term into the Lagrangian which controls the quantum and thermal contribution 
to the effective action. For the bosonic modes, we insert the bosonic mass-like 
term R/ : b(p)4 > (p) ( I ) (~p) m to the Lagrangian where RkB{p) is an arbitrary function 
satisfying the following condition: 

RkBip) ~ k 2 for p 2 < k 2 , 
R kB (p) ~ for p 2 » k 2 . (4) 



At the infrared scale (k = 0), Rk is zero for any momentum, and becomes the 
quantum effective action. At the ultraviolet scale (k = A), all modes acquire large 
masses, m 2 — > m 2 + aA 2 where a is a proportional coefficient in Rkip), Ra{p) — ctA 2 . 
As a result, R^ prevents the propagation of all the fluctuation, and lyi becomes the 
classical action. Thus one finds that the conditions for the 1^ are satisfied. At an 
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intermediate scale, Rk(p) suppresses propagation of the modes whose momentum are 
smaller than k. Thus we can interpret the r k as an effective action which contains the 
fluctuations whose momentum scales are larger than k. The above procedure can be 
extended to the fermionic modes by introducing the mass-like term ip(p)RkF(p)' l P(p) 
where R^f has Dirac indices. (See, e.g., RefPO 1 ). 

In the present work, we choose the so-called optimized cutoff function^ whose 
form is 

R kB (p) = (k 2 -p 2 )9{k 2 -p 2 ) for Bosons 

( fk 2 \ -, ( 5 ) 

Rwip) = f\ \ ~^Tn — 1 6(k 2 — P 2 ) fo r Fermions . 
\\p J 

The evolution of the with scale k is described by the FRG equationPS 



k^ - -Tr 
k dk ~ ir 



dkRkF 



R kF + r^ 



+ -Tr 

2 

F 



dkRkB 



RkB + /f 2 > 



(6) 



where the trace runs over Dirac, color and flavor indices as well as momentum. The 
superscript (a, b) in rjf"'^ denotes a(b) times fermionic (bosonic) functional derivative 
of the Jfc. The first and second terms of r.h.s of © are the fermionic and bosonic 
contributions respectively. 

To solve the FRG equation of the scale dependent effective action ©, we need 
scale (k) dependencies of the 2-point functions rj?' ^ and r^' 2 \ In general, the 
FRG equation for the n-point function depends on up to the (n + 2)-point function 
and the set of the FRG equations for n-point functions forms an infinite tower of 
the equations. In practice, we introduce some approximations in order to close the 
tower of equations up to finite n-point functions. In the present work, we apply the 
so-called local potential approximation (LPA) to the effective action. We assume the 
form of the scale dependent effective action as 

F k LPA = - g s (a + ij 5 f ■ vr) - g d ifr/o + "7o]V> + \{d^) 2 + ^{d^) 2 + ^(d^) 2 
- U k (a 2 + TT^,ip) - ca , 

(7) 

where U k is a scale dependent effective potential such that it agrees with a usual 
effective potential at k = 0. Except for the small explicit breaking term ca, which can 
be treated as a perturbation for the symmetric effective action, the Lagrangian ([2]) 
respects the chiral symmetry. To a good approximation, we can regard U k as a 
function of a 2 + n 2 and (p. 

Here we neglect the scale dependence of the wave function renormalization for 
simplicity. It would be important for the discussion of the possibility of the inhomo- 
geneous condensation phase but it is out of the scope of this work. 

We also neglect the evolutions of the Yukawa coupling g s and the density coupling 
gd- The inclusion of the scale dependencies of these couplings might change the global 



G 



structure of the phase diagram such as the position of the QCD critical point or the 
slope of the first-order phase boundary of the chiral phase transition. Near the QCD 
critical point, however, the contributions from the light bosons are enhanced and the 
scale dependence of g s and would not change the local phase structure, i.e., the 
shape or size of the critical region. Thus the scale dependence could be neglected 
when the critical region is only concerned. 

Now the 2-point functions in the r.h.s of ([6]) is given by the functional derivatives 
of the LPA effective action: 



(2,0) 



r, 



(0,2) 
k a,b 



dtpdip 

d 2pLPA 



(8) 



Using LPA, we can obtain an evolution equation for the scale dependent effective 
potential. Substituting the LPA ansatz ([7]), ([8]) and the optimized cutoff function 
© into ©, we obtain the flow equation for the effective potential: 



dU k 



k l 



dk 12vr 2 







- 2N f N c 


— tanh ( 
_E q \ 



'E q + (n + g d ip)\ 1 



2T 



+ 



E„ 



tanh 



E q -{n + gd¥>) 



2T 



(9) 



where the scale dependent particle energies E a (a = q, ir and ±) are defined as 



E 2 a = k 2 + Ml 

M 2 q =g 2 a\ M 2 = —, 
a 



M 



(10) 



U" + U ± J(U" -U) 2 + 4U' 2 



i=± 



Here we have used the following notation 



aa dip 



(11) 



M 2 and M 2 are understood as the squared curvature masses of the ir and the linear 
combination of a and (p. 

Since the scale dependent effective action at the UV scale agrees with the classical 
action © , the initial condition of the flow equation Q reads 



U A (a 2 + vf 2 , <p) = a(a 2 + vf 2 ) + b(a 2 + vr 2 ) 2 + 



(12) 
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a /A 2 


b 


c/A 6 


M^l A' 


A MeV 


9s 


0.3224 


0.25 


0.0045 


0.1 


900 


3.2 



Table I. Table of the initial parameters for the numerical calculations 



By solving the flow equation from k = A to k = 0, we obtain the effective potential 
with the fluctuation effects incorporated. The expectation values of the a and (p are 
determined by minimizing the effective potential at k = 0: 



dU 
da 



dU 

dip 



. 



(13) 



We obtain the effective potential at various temperatures and chemical potentials, 
and the phase structure is determined by the chiral order parameter behavior. 



§3. Numerical results 



In this section, we show numerical results. One of our concerns is the sensitivity 
of the phase structure to the strength of the mixing between the chiral condensate 
and the quark- number density, i.e., the magnitude of the density coupling g^. 

First we calculate the phase diagrams with g& being varied. We show the position 
of the QCD critical point slightly moves to a higher temperature direction with 
increasing g^. 

Next we describe the critical region which is associated with the QCD critical 
point. Finally we calculate the curvatures of the effective potential. We show that a 
kind of level repulsion between the curvature masses occurs around the QCD critical 
point at finite g^. This level repulsion causes an expansion of the critical region. 

3.1. Calculation method and parameter fixing 

We solve the flow equation (|UJ) numerically on a grid in the field space. Since 
both a and (p can take non-vanishing expectation values, we have to prepare the 
Uk values on the two-dimensional grid of discrete field values in the a and <p space. 
Derivatives of the effective potential used in the numerical calculation are given by 
finite differences. This reduces the flow equation to a coupled ordinary differential 
equations which can be integrated out with a standard method. We solve it by using 
the fourth-order Runge-Kutta method. 

There are few works using a 2-dimensional grid method for dynamical fields 
and, only very recently, it has been developed by several researchers including one 
of the present authors on the studies of QCD phase diagram with isospin chemical 
potentiaPD'ES 1 as we \\ as th e phase diagram of 2-color QCDP^ 1 In comparison with 
the Taylor method in which the effective potential is locally evaluated around a scale 
dependent minimum, the grid method captures the global structure of the effective 
potential in the field space and is well suited to investigate the critical region with 
which the first order phase boundary is connected. 

In the present calculation, we choose the ultraviolet scale as A = 900MeV. Our 
effective model has six parameters a,b,c, g s , M v and g^. The first four parameters 
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have been determined to reproduce the vacuum physical values such as the pion 
mass, sigma mass, pion decay constant and constituent quark massF^ We can not, 
however, determine the remaining two parameters M v and g^ because at present 
any experimental data for the collective excitation corresponding to the particle- 
hole mode (phonon mode) in quark matter is not available. In principle we could 
determine these parameters from microscopic theories such as QCD itself or the NJL 
model. However the main purpose of this work is to examine the sensitivity of the 
critical region to the mixing between the a and then we treat and as free 
parameters and examine responses to their variation. Here we fix = 0.1/1 2 and 
change g^. The phase structure is determined by the relative strength g v /M^. For 
example, the ordinary results of the QM model are reproduced by taking g v /A'l^ — > 0. 
Our choice of the parameter set is summarized in Table HI 

3.2. Phase structure 

Before studying the critical region, we start with an analysis of the phase struc- 
ture with varying gj. In Fig. [TJ the phase diagrams are shown for four fixed values of 



100 r 
80 - 
> 60 

S 

H 40 

20 

- 



> 



100 p 
80 
60 
40 

20 ■ 

- 



g d = o 



z 



> 



180 210 240 270 300 330 360 
|x[MeV] 



g d = 0.2g s 



1 



180 210 240 270 300 330 360 
|x[MeV] 



100 
80 
60 
40 
20 





1 

100 
80 
60 
40 
20 




g d = o.ig s 



z 



180 210 240 270 300 330 360 
ufMeV] 



g d = -3g s 



180 210 240 270 300 330 360 
|x[MeV] 



Fig. 1. Phase diagram for four different gd (0, 0.1g 3 , 0.2g s and 0.3g s )- The solid line denotes the 
first-order phase boundary while the blue closed circle represents the critical point. A green 
line and black circle in g d = 0.3g s denote second first-order phase boundary and its critical end 
point. 



gd- In the present work, we have generalized the previous results in the zero density- 
coupling case^2J to the finite couplings. We find the QCD critical point for each 
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value. The locations of the QCD critical point are, in unit of MeV, at (T cp ,fj, cp ) ~ 
(15, 303), (16, 303), (28, 300) and (42,284) for g d = 0,0.1^,0.2^ and 0.3ft,, respec- 
tively. For the temperature below the critical point, we find the first-order phase 
boundary which separates a chiral restored phase (higher fi) and a broken phase 
(lower fi). At the critical point, the phase transition is second order. Beyond the 
critical point, the phase transition becomes crossover and the pseudo critical line is 
not shown in the phase diagram. 

The position of the critical point slightly rises to a higher temperature direction 
with increasing gd- This is an opposite behavior to the well known results in the 
NJL model with the vector interactionP3fi3JIl}SZl | This is because our effective 
model corresponds to the attractive vector interaction case (Gy < 0). In fact, a 
similar behavior is found in the mean field approximation of the NJL model with an 
artificial attractive vector interactionPS 1 

In the figure for <jy = 0.3g s , one sees another critical point at (T cp ,{i cp ) ~ (80,351) 
MeV from which another first-order phase boundary is extended to the lower-T di- 
rection. This should be an artifact of our model. This additional critical point and 
its associated first-order phase boundary appear only with large couplings such as 
gd = 0.3g s and they belong to a different branch from the QCD critical point. The 
magnitude of the chiral condensate is sufficiently smaller than the vacuum expecta- 
tion value and jumps from about 10 MeV to a few MeV across the first-order phase 
boundary. At g^ = 0.3g s , the attractive force between two quarks is too strong and 
causes another stable minimum of the effective potential in the high density and 
small chiral condensation region. 

In Fig. [21 the effective potential contour is shown for several values at temper- 
ature and chemical potential close to the critical point. Since the phase transition 
is second order at the QCD critical point, there is a flat direction on the effective 
potential around the minimum. 

For gd = 0, ip is decoupled from the chiral modes then the mixing term such as 
a 2 <p does not arise in the effective potential. As a result, we can see the flat direction 
is parallel to the a axis. 

On the other hand, for gd ^ 0, there is a mixing between a and ip and the 
expectation value of ip is not zero. 

The flat direction is tilted to the direction of the linear combination of a and ip 
directions and the a direction is no longer special. These behaviors are also found 
in a calculation on the NJL modelPD The angle between the flat direction and the 
a direction becomes larger with increasing gd as the mixing between the a and ip 
becomes stronger. 

3.3. Susceptibility 

The information on the size of the critical region is important for experiments 
in search of the QCD critical point. We define the critical region by using the quark- 
number susceptibility \ q given as the response of the quark number density n q to 
ft: 



dn q 



djj? 



(14) 
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40 60 80 100 40 60 80 100 



o [MeV] [MeV] 

Fig. 2. Contours of the effective potential in the a and ip space near the critical point for varying 
gd- The solid points show the minimum of the effective potential. 



where the pressure P is given by P(T,fj,) = — U{(Jq,^)q). For the comparison of 
the shape and size of the critical region for different g^, it is convenient to use a 
normalized susceptibility Xqi i- e -> the ratio to that of the mass-less free quark gas 

v free . 
Xq 

free = ™ c Nf 

Xq ~ 6 

Xq Xq/ Xq 

In Fig. [3l Xq a t 9d = 0.2<7 S is shown as a function of quark-chemical potential for 
three different temperatures. For temperatures below the critical temperature T cp , 
there is a discontinuity in Xq a t the first-order phase boundary. For temperatures 
above T cp , Xq is continuous as a function of fj, and has a peak at the crossover phase 
boundary. The peak height b ecomes higher with T approaching to T C p. At T — ^cp; 
the peak height of Xq diverges. Near the QCD critical point, we find the same 
behavior of the Xq f° r an 9d- 

In Fig. HI the contour plots of Xq with two fixed values, Xq = 3, 6, are shown 
for varying g^. The criticality near the QCD critical point is characterized by the 
enhancement of Xq- I n t ne following, we define the critical region by Xq > 3. 

For g<i = 0, we have the critical region which has 10 MeV extent both in the 
T and \x directions. The size of the critical region for = is consistent with the 
previous worlPO 1 anc j the width of the critical region is shrunken in the [i direction 
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Fig. 3. Plots of the normalized quark-number susceptibility \q — Xq/Xq C ° with chemical potential 
at Qd = 0.2g s . The red/blue line represent the \q a ^ T below/beyond the T cp . The green line is 

at 1 = Te n . 
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Fig. 4. Contours of the normalized quark- number susceptibility x!J = Xq/Xq CC around the critical 
point for \g = 3 and 6. The red line denote the first-order phase boundary. The blue solid 
circle represent the QCD critical point. 



by almost one order of magnitude compared with the mean-field calculation.^ 

For finite gj, the critical region is deformed from that at g& = 0. We observe 
the expansion of the critical region to the crossover phase transition direction with 
increasing g^. The shape and size of the critical region are quite sensitive to g^ such 
that it shows drastic deformation while the position of the QCD critical point moves 
only slightly. At = 0.3g s , the critical region is elongated to the cross over side 
twice as large as that for = 0. 
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This expansion can be partly understood by looking at the curvature masses of 
the linear combination of the a and ip, as demonstrated in the following. 

3.4. curvature masses 

In order to give the understanding of the expanding behavior of the critical 
region, we calculate the curvature masses of the linear combination of the a and 
ip. The masses M + and M_ are defined as a bigger and smaller eigenvalue of the 
curvature matrix, respectively: 



M 



( d 2 u 

I (9(7(9(7 

I d 2 U 

\dipda 



d 2 U \ 

d 2 U I 

dtpdtp / 



(16) 



In Fig. the curvature masses M± at T 
chemical potential for = (left) and g& - 



- T cp are shown as functions of quark- 
0.2g s (right). 
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Fig. 5. Plots of the curvature masses M± near the critical point with chemical potential. The left 
and right panel are at g& = and ga — 0.2g a respectively. 

At gd = 0, there is no mixing between the a and (p then the non-diagonal parts 
of the mass matrix (|16p are always zero, and thus we can identify the a and cp modes 
individually. Since the <p is decoupled from the matter sector, the curvature mass of 
the p is constant (~ 400 MeV) with chemical potential varied while the sigma mass 
decreases below the mass of the (p and vanishes at the critical point. Beyond the 
critical point, the a mass increases monotonically. 

In contrast, the non-diagonal parts of the mass matrix (|16p do not vanish for 
gd so that the a and <p are mixed. At g^ = 0.2g s , M_ stays constant (~ 400 
MeV) while M + decreases below the critical chemical potential. When M + gets 
close to M_, M_ starts to decrease and vanishes at the critical point. We observe a 
kind of level repulsion near the critical point between M + and M_ which is familiar 
in a perturbation theory in quantum mechanics!® Beyond the critical point, M_ 
increases until 340 MeV and becomes almost constant. M + increases again while 
M_ increases and approaches a constant value. We found these behaviors of the M + 
and M_ generally for ^ 0. 

Note that the softening behavior shown here does not correspond to the dynam- 
ical softening of these modes. In order to see the dynamical softening of the modes, 
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one has to calculate spectral functions of the a or (p modes in which one would find 
a dissipative zero mode in the space-like region of the spectral function 
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Fig. 6. Contours of the smaller curvature mass M- for varying gd- The contours are shown for 
M- = 200 and 300 MeV. The red lines denote the first-order phase boundary. The blue solid 
circles represent the QCD critical point. 



In Fig[6l we show a contour map for M_ with varying g^. For each g^, M_ is close 
to zero in the vicinity of the critical point. The softening of the linear combination 
of the a and ip causes the critical behavior on the susceptibility. 

In the following, we call the area in which M_ < 200 MeV the softening region. 
For each gd, we can see the softening region almost coincide with the critical region, 
i.e., the region in which Xq > 3 in Fig. UJ The softening region is also expanded to 
the crossover phase transition direction with being increased. 

The level repulsion of the curvature masses precisely corresponds to the expand- 
ing behavior of the critical region. The curvature mass M_ is reduced when the 
density coupling or the non-diagonal parts of the mass matrix ()16p is finite. The 
non-diagonal parts of (|16j) become larger as gd increases. As a result, the level re- 
pulsion between M + and M__ becomes larger and the region with small M_ also 
gets expanded. It simply follows that the inclusion of the mixing between the a 
and tp and does not depend on the details of interactions, for example, repulsive or 
attractive force. 
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§4. Summary 

We have examined the influence of the mixing between the chiral condensation 
and the quark-number density on the critical region around the QCD critical point. 
The fluctuations of the quark-number density as well as the chiral modes were taken 
into account through the functional renormalization group (FRG) equation. 

We have extended the quark meson model to a new effective model which is ap- 
propriate for the correct description of the dynamics near the QCD critical point. In 
addition to the quarks and chiral modes (cx and tt) , the new effective model contains 
a new field ip which corresponds to the quark-number density. The magnitude of the 
mixing between the a and ip was determined by a density couping g^. 

We derived the flow equation for the scale dependent effective action with the 
fluctuation of the density field at a finite temperature and quark-chemical potential. 
The flow equation for the scale dependent effective potential was solved numerically 
on a two-dimensional grid of the chiral condensate (a) and quark-number density 

(¥>)■ 

We found that the QCD critical point exists for varying ga while its position 
slightly moves to the higher temperature and lower quark-chemical potential direc- 
tion as gd increases. The quark-number susceptibility was enhanced around the 
critical point. 

In contrast, we have shown the shape and size of the critical region was quite 
sensitive to the magnitude of the mixing term. We have shown that the critical 
region is enlarged to a lower chemical potential and higher temperature direction 
when we increase the density coupling g^. 

In order to confirm the expanding behavior, we also calculated the curvature 
masses of the linear combination of the a and tp. It turned out that the critical region 
was almost identical to the softening region where the curvature of the effective 
potential becomes small. The softening of the linear combination of the a and p 
causes the critical behavior. 

We have shown that the origin the expanding behavior of the critical region and 
the softening region can be partly attributed to the level repulsion of the curvature 
masses given by the linear combination of the chiral condensate and the quark- 
number density. The large level repulsion at large ga makes the critical and softening 
region wider. The level repulsion is a generic consequence of the mixing between the 
chiral condensate and the quark-number density. 

As a future work, it is interesting to investigate higher moments of the quark- 
number density or the charge-number density including the quark-density coupling 
and the fluctuation of the quark-number density. The higher moments are considered 
to show the stronger criticality^ 1 and would be more sensitive to the mixing than the 
quark-number susceptibility. In the present work, we have considered the simplest 
case of the coupling and treated its strength as a free parameter but it should be 
determined by a more microscopic model or experiments. The exact determinations 
of the strength or detailed form of the mixing also awaits future studies. 
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